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We studied magnetic properties of the double exchange (DE) model with S — 1/2 localized 
spins at T = 0, using exact diagonalization in the framework of the dynamical mean field theory. 
Obtained phase diagram contains ferromagnetic, antiferromagnetic and paramagnetic phases. 
Comparing the phase diagram with that of the DE model with classical localized spins, we 
found that the quantum fluctuations of localized spins partly destabilize the ferromagnetism and 
expand the paramagnetic phase region. We found that phase separations occur between the 
antiferromagnetic and paramagnetic phases as well as the paramagnetic and ferromagnetic ones. 
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§1. Introduction 

Recently perovskite manganese oxides doped with 
divalent alkaline metals, such as Lai_a;Sra;Mn03 and 
Lai_2;Caa;Mn03, are studied extensively. In these ma- 
terials the decrease of 3d electrons caused by the doping 
with divalent metals leads to drastic changes of mag- 
netic and transport properties. For example, the mother 
material LaMnOa is insulating and has an antiferromag- 
netic (AF) order at low temperatures. This AF insula- 
tor is transformed to a ferromagnetic (F) one by doping, 
and the F state becomes metallic for x > 0.20.^-' The 
orthorhombic crystal structure changes into a rhombo- 
hedral one with higher symmetry in the metallic phase. 
The so-called colossal negative magneto-resistance shown 
by these materials is attracting intense interests due to 
its potential applications.^^ 

Zener proposed the double exchange (DE) mechanism, 
which describes essential physics in these materials. In 
this paper we study the simplest model that contains this 
mechanism, i.e. the DE model. In this model three t2g 
electrons on a Mn ion are treated as a localized spin and 
the two-fold degeneracy of eg orbitals is neglected. The 
Hamiltonian is defined as 

<i,j>,a i 

Here t denotes the hopping amplitude between the near- 
est neighbor sites, Cia (c^) the annihilation (creation) 
operator of an eg electron at the site i with spin a, and 
Si and Si are the localized spin and the spin of the Cg 
electron at the site i, respectively. Itinerant eg electrons 
are assumed to occupy a single orbital and interact with 
the localized spins through the Hund coupling Jjj. 



Anderson and Hasegawa^-* showed that in the limit of 
strong Hund coupling the hopping amplitude of electrons 
between the sites i and j is reduced as Uj = tcos(0y/2) 
when the angle between the localized spins Si and Sj 
is 9ij. The hopping amplitude is largest when 9ij = 0, 
which implies that the F state minimizes the kinetic en- 
ergy. This is the DE mechanism for the ferromagnetism 
proposed by Zener. The competition between the DE 
mechanism and the antiferromagnetic super-exchange in- 
teractions between the localized spins was studied by de 
Gennes.^^ Kubo and Ohata studied the DE model in the 
strong coupling limit. 

Recently Furukawa studied this model with classical 
localized spins using the dynamical mean field theory 
(DMFT), which is exact in infinite dimensions. In 
DMFT the effective action for the system is equivalent 
to that for a single impurity problem, and the momen- 
tum dependence is absent in the self-energy of the one- 
particle Green function. These features of the DMFT 
make the calculations of various properties tractable. He 
obtained the thermodynamic as well as the transport 
properties of the model employing a semicircular as well 
as a Lorentzian density of states for non-interacting elec- 
trons. The DE model with classical spins was studied 
also in finite dimensions with numerical methods. Ow- 
ing to these studies the properties of the DE model with 
classical localized spins are fairly well understood. The 
DE model with quantum localized spins is also studied. 
The spin wave spectrum^^^ and the ground state phase 
diagram"'^^' ^''^ in one and two dimensions were studied by 
numerical diagonalization. Variational arguments^^' 
show that the fully polarized state is destabilized in a 
certain density region even in the strong coupling limit. 
Electronic states were also studied by a many-body CPA 
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approximation.^^) Due to the difficulty of the strong 

correlation effects, wc have still much to understand the 
effects of quantum fluctuations in the DE model. 

In this paper, we study the DE model in the cases of 
= 1/2 and the classical (the limit 5 = oo with finite 
JjjS') localized spins in infinite dimensions. Comparing 
the magnetic phase diagrams in both cases, wc clarify 
the influence of the quantum fluctuations of the local- 
ized spins on the magnetic ordering. We use the exact 
diagonalization method in the framework of the DMFT 
in order to study the case of S* = 1/2. Wc employ the 
Gaussian density of states for non-interacting electrons, 
which is exact for the simple hypercubic lattice in infinite 
dimensions, 



^o(e) 



1 



(1.2) 



The contents of this paper are as follows: In section 
2 we summarize the method of studying the DE model 
in infinite dimensions. In section 3 we report the phase 
diagram and the magnetic susceptibility for the case of 
5* = oo. In section 4 we show the results for the case of 
5=1/2 and compare them with those for 5 = oo. In 
section 5 we summarize the results and give perspectives 
for future study. 

§2. Method 

We employed the DMFT to study the DE model (1) in 
infinite dimensions. Though the method is now a stan- 
dard one,^*) we briefly describe it in this section for read- 
ers who are not familiar with the method. In infinite 
dimensions, the infiuence from neighboring sites can be 
replaced with a dynamical mean field (Weiss field). The 
problem is hence reduced to a single impurity problem 
in the Weiss field. The self-energy of the Green function 
does not depend on the momentum. ^^'^"^^ The Green 
function is written as 

G,0^,iiJn) = -. ^ W-^' (2-1) 

where Scr(iw„) is the self-energy and a denotes the spin 
suffix. It depends on the momentum k only through ej^. 
In the following text we omit the variables k, iuj and cr 
in G, S etc. for simplicity. The effective action of the 
DE model at a temperature T = (3~^ is 

„ Jo Jo 



Jh r drS{T)-s{T), 
Jo 



(2.2) 



where 'if'^(T) and ^'^(t) are Grassmann variables rep- 
resenting the electrons at the impurity site, S(r) the 
impurity localized spin, and Gq^ is the Weiss function 
that represents a time-dependent mean field. This Weiss 
function is related to the Green function G and the self- 
energy S as = G"^ - S. 



We obtain the Weiss function Qo and the Green func- 
tion G self-consistently through the following procedures. 

A. Given Gq^, then calculate G employing the effective 
action. 

B. Calculate S from G through the relation Qq^ = 
G^i - S. 

C. Calculate G from E as 

^o(e) 



G{iuj„ 



de- 



(2.3) 



iuJn — e — /i — S 

where Dq{€) is the density of states given by (1.2). 
D. Obtain Gq^ again by using Gq^ = G~^ — S. 
Wc repeat the procedures from A to D until a self- 
consistent solution is obtained. 

The main difficulty in these procedures lies in the pro- 
cedure A, where we must solve a quantum many-body 
problem. For finite S, we cannot accomplish this proce- 
dure exactly and hence we need an approximation. The 
Weiss function is expressed as 



= iu!„ + fi + 'y^ 



duj' 



p=2 ^P'^ 



(2.4) 



where /^ai'jj) = Y^'^ V^cr^{uj — ipa)- In this study we 
approximate Goaiii^n) with Goai^^n) which is the Weiss 
function for a finite number of discrete orbitals; 



Ns y2 



p=2 



(2.5) 



The parameters Vp„ and ip„ are determined self- 
consistently as to minimize 



, nmax 

R= —r E \Go{i<^n)-^ - Gr^\ioJn)-'\, (2.6) 

n—0 

where rimax labels the maximum frequency employed in 
the calculation. When the Weiss function is approxi- 
mated as above, the problem is identical to solve the 
following impurity Anderson Hamiltonian, 

iVg JVg 

W = E epcra+^ap^+ ^ Vp^{a+^d„ + d+Up^) 

p>2,<j P>2,cr 



+edY^d+da - JrS • s, 



(2.7) 



where a^^ [dpa) and d^ {da) are the creation (annihila- 
tion) operators of the orbital p and the impurity site, re- 
spectively, and s" = ^ Ylaa' dt'^a cr'da'- Wc obtain the 
ground state of this Hamiltonian with use of the exact 
diagonalization method, and calculate the Green func- 
tion. The kinetic and interaction energies in the ground 



Magnetic Order in the Double Exchange Model in Infinite Dimensions 



3 



state |G.S.) are obtained with use of the expressions,^^^ 
E^in = (G.S.I ^ Vp{a+,d^ + dtap,) |G.S.) (2.8) 

P>2,<T 

and 

^int = (G.S.I -JhS-s|G.S.). (2.9) 

Although the size of a localized spin is 3/2 in actual 
manganese oxides, we adopt the; case of 5 = 1/2 in order 
to emphasize the quantum effects as well as for simplicity. 
The approximation improves with increase of the number 
A^g — 1 of the assumed orbitals. The present study was 
restricted to Ng < 8 due to technical reasons. 

We also calculate the case of 5 = oo for the sake of 
comparison. In this case the procedure A is easily ac- 
complished becaiisc the effective action is that of non- 
interacting electrons. The partition function in the case 
of 5 = oo is given by 

Z = Jdn^J Ylv-^lV^^e^pi-SiS)) (2.10) 

with 

S{S) =Y.f f dT2K{ri)Qoa{ri - r2)*.(T2) 
" Jo Jo 

- Jh / drS ■ s(r), (2.11) 

Jo 

where S is a classical vector of size S and fig denotes 
the volume angle of the direction of the localized spin. 
The Green function is given by 

Ga(iw„) = ^ y df2sP(S) {Qoa{i^n) + JhS • s(ia;„))"' , 

(2.12) 

where P(S) is the Boltzmann weight for the configura- 
tion of the localized spin, 

P(S) = / n 'DKV-^, exp (-5(S)) . (2.13) 

§3. Results for S = oo 

In this section we report the temperature dependence 
of the magnetic susceptibility and the phase diagram at 
a very low temperature in the case of S' = oo. 

In infinite dimensions the susceptibility depends on 
the wave vector q only through A(q),^°) where X(q) is 
defined by 

1 

A(q) = -^cosft. (3.1) 

i=l 

Here d denotes dimensionahty of the space. The F 
and AF susceptibilities correspond to the cases with 
X{q} = 1.0 and —1.0, respectively, and a case where 
— 1.0 < A'(q) < 1.0 corresponds to an incommensu- 
rate order. Figure 1 shows the magnetic susceptibility 



X for JhS" = 5.0 at the electron density n = 1.0. Figure 

1(a) exhibits the wave- vector dependence of the suscep- 
tibility at various temperatures. The susceptibility has 
maximum at X{q) = —1.0 at all temperatures and it 
increases with decrease of temperature. In Fig. 1(b) the 
inverse of the AF susceptibility is shown as a function of 
temperature. The susceptibility obeys the Curie- Weiss 
law and diverges at a finite temperature . As a result 
the AF order is expected to appear at temperatures be- 
low Tn = 0.044 for JuS = 5.0 and n = 1.0. Next, the 
susceptibility for J^S = 5.0 and n = 0.8 is shown in Fig. 
2. In this case the F susceptibility is largest, and obeys 
the Curie- Weiss law. The F order is expected to appear 
at temperatures below Tq = 0.062. 

We obtained a phase diagram at a very low tempera- 
ture (/? = 400) assuming the uniform or a two sublatticc 
structure. The result is shown in Fig. 3, where the AF, 
F and paramagnetic (P) phases arc present. For ri = 1 
the AF state is stable for any positive Jjj. The AF phase 
spreads out to a small region close to the line n = 1 for 
weak Hund coupling. 

The F phase occupies a large region where n <, 0.5 
and JhS ^ 1. This F state has almost fully polar- 
ized magnetization. In Fig. 4, we show the spontaneous 
magnetization for n — 0.4 at /3 = 400 as a function 
of JylS. The magnetization starts increasing continu- 
ously at Jh5' = 0.9 and rapidly approaches the limiting 
value of 0.98mo where mo is the value of the full polar- 
ization. In this case the transition to the F state is of 
the second order. The P phase appears for JhS ^ 3 
and it spreads to a wide density region for weak cou- 
pling. Phase separation (PS) occurs between the AF 
and P phases for Jh^* ^ 3 and also between the AF and 
F phases for J^S ^ 3. General features of the phase 
diagram is similar to those obtained previously in infi- 
nite dimensions'^' and roughly understood as a result 
of the competition between the AF coupling, which is 
a second order effect of the hopping, and the ferromag- 
netic DE mechanism. Since the AF coupling decreases 
with increasing J^S as t^/JiiS, the relative strength of 
the F coupling, which is proportional to t, increases with 
Jh. Occurrence of the PS for small doping (n ~ 1) may 
be understood as the instability of a canted AF state. 

Since we assumed a two-sublattice structure there re- 
mains the possibility that we missed magnetic phases 
with more complicated spatial structures, e.g. those 
with incommensurate (IC) wave vectors. Such IC order 
was predicted previously to occur for intermediate dop- 
jj-^g 22. 24) r^Yic obtained P phase might be an IC phase in 
reality. In fact we found the cases where x corresponding 
to an IC order diverges at higher temperatures. Figure 
4 shows X for JjiS = 1.0, n = 0.66 and P = 110. We 
observe a sharp maximum of x at X{q) ~ —0.9 and it 
diverges at a critical temperature. It is certain that an 
IC phase exists below the critical temperature. Study of 
the phase diagram which contains IC phases is left for a 
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future study. 

§4. Results for 5= 1/2 

Quantum Monte Carlo method is usually used in the 
framework of the DMFT for studying the thermody- 
namic properties at finite temperatures. This method, 
however, suffers from a serious negative-sign problem 
in the study of the DE model with quantum localized 
spins. In this section, we therefore studied only the 
ground state properties for the model with S = 1/2 lo- 
calized spins, employing the exact diagonalization tech- 
nique. We report the results and compare them with 
those for the case of 5* = oc'. 

Figures 6(a) and 6(b) show the groimd state phase di- 
agram obtained with Ng = 6 and 8, respectively. We 
should extrapolate the results of finite A^s to the limit 
A^S = oo in order to obtain the exact result. We observe 
that the two phase diagrams are almost same except for 
the boundary of the F phase. Difference between results 
for Ns = 6 and 8 is still large on the boundary of the F 
phase and hence the obtained boundary may not repre- 
sent that of A^'s = oo correctly. Nevertheless these results 
exhibit contrast to that for the case of S' = oo in several 
points. 

The AF phase appears in a region close to the line 
n = 1 and its boundaries for A^g = 6 and 8 agree very 
well. It is remarkable that the AF phase is much more 
stabilized for S = 1/2 than for S' = oo in the weak 
coupling region. It extends to n ~ 0.7 and, as a result, 
the region of PS narrows for small Jjj in contrast to the 
result for S' = oo. The present result shows that the 
AF phase is stable for 0.8 <, n < 1 even for Jjj = 0.1, 
the smallest Jjj studied, though the system should be 
paramagnetic for all < n < 1 for Jjj = 0. 

The F phase appears for J^S ^ 4 in a wide range 
of the density, i.e., n 0.7. Results for A^g = 6 and 
8 agree well in this region. On the other hand, the F 
phase is strongly reduced for JhS* ^ 3 by increasing A^g 
and almost disappears for J-^S <; 1 in Fig. 6(b). In 
comparison with the S = oo case, the F phase region for 
A^g = 8 is fairly reduced. We may need further study 
with larger Ns to confirm this strong reduction of the F 
phase. 

Between the AF and F phases we obtain a P phase in 
a wide density range, and PS occurs between the AF and 
F phases as well as between the P and F phases. The 
P phase persists to large Jh, i.e. J^S ~ 5 in contrast 
to the case of 5 = oo where the P phase terminates at 
JfjS' — 3. We should note that the obtained P phase 
might be replaced by phases with more complicated spa- 
tial structures, e.g. IC phases, in more comprehensive 
calculations, which include larger-sublattice structures. 

Figure 6 shows the electron density n as a function 
of the chemical potential /x for J^S = 5.0 and A^s = 8. 
We see jumps of n at two values of /x, which indicate 
the occurrence of PS's. The density region of the PS 



between the AF and P phases narrows with decreasing 
Jh in contrast to the case of S' = oo. This is caused 
by the expansion of the AF phase in the weak coupling 
region. The PS occurs in a very narrow density region 
between the P and F phases for Jh-? ^ 2. This PS was 
not found in the case of = oo. 

In Figs. 7(a) and 7(b), we show the kinetic energy 
and the magnetization of conduction electrons, respec- 
tively, as functions of Jh^ for n = 0.4. The kinetic 
energy [Fig. 8(a)] increases with Jjj very rapidly in the 
P phase, while it shows a saturation in the F phase for 
JrS > 3.5. This behavior is a characteristic of the ferro- 
magnetism caused by the DE mechanism and indicates 
that the scattering by the Hund coupling ceases to be 
operative by alignment of the localized spins. ^"-'^ In Fig. 
8(b), the F phase starts with a finite magnetization at 
JllS = 3.5, which indicates a first-order transition from 
the P phase. The occurrence of PS between the F and P 
phases is related to the first-order character of the phase 
transition. The magnetization increases with Jjj and 
apparently approaches to a finite value nearly equal to 
0.6mo, where toq denotes the magnetization of a fully po- 
larized state. It is noticeable that the saturated moment 
is relatively small in comparison with the result for the 
15 = 00 localized spins. This result gives evidence that 
the quantum fluctuations of the localized spins destabi- 
lize ferromagnetism. Though the fully polarized state 
cannot be realized for the Gaussian density of states, the 
reduction of the magnetization mostly comes from the 
quantum fluctuations of the localized spins. 

We also studied the electronic states in the groimd 
states. The local density of states (DOS) is calculated 
from the Green function obtained in the procedure A 
in section 2. Because we treat the system with a finite 
number of discrete orbitals, the obtained DOS is com- 
posed of finite number of sharp peaks. Figures 9(a-d) 
exhibit examples of the DOS obtained by calculations 
with A^s = 8, where an artificial imaginary part O.li is 
added to S(w) for the purpose of illustration. Figure 
9(a) exhibits the local DOS in an AF ground state for 
JuS = 5.0 at n = 1.0. We can observe characteristic fea- 
tures of the DOS in the AF ground state in the strong 
coupling region (in spite of the discretized structures due 
to the approximation). In this case DOS is composed of 
two sub-bands each of which corresponds to electronic 
states with parallel and antiparallel spins with localized 
spins. Two sub-bands are separated from the Fermi level 
by a broad energy gap, which indicates that the ground 
state is insulating. The lower DOS for up spin has a large 
weight compared to that for down spin, since spin is al- 
most fully polarized upwards on this site. Figure 9(b) 
shows the local DOS of a F ground state for J^^S = 5.0 
at n = 0.78. Also in this case the DOS is composed 
of two sub-bands separated by an energy gap (the main 
part of the upper sub-band is out of range of the fig- 
ure). The Fermi level is inside of the lower sub-bands. 
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which indicates the metaUic partially polarized ground 

state. The DOS for down spin is much smaller than 
that for up spin, which implies that the spins are po- 
larized upwards. The width of the lower sub-bands for 
both spins are almost same, which reflects the mixing 
due to exchange scattering. The DOS in the AF ground 
state for J^S — 0.5 at n = 1.0 is exhibited in Fig. 9(c). 
In this case the DOS is apparently composed of a single 
band centered around the Fermi level. It is expected that 
there is an energy gap caused by the AF order and/or 
Umklapp scatterings. However it is exhibited as a dip 
at the Fermi level. This is only due to the width added 
artificially. The DOS for down spin below Fermi level 
is smaller than, but comparable to that for up spin due 
to small sublattice magnetization. We show in Fig. 9(d) 
the DOS for J^S = 0.5 and n = 0.73 where the ground 
state is antiferromagnetic. In this case the dip located 
at e = /i + 0.3 may correspond to the gap due to the AF 
ordering. The Fermi level lies very close to the top of a 
peak. These results suggest that the ground state is a 
metallic state with AF long range order. 

§5. Summary and future work 

We studied the DE model with S = 1/2 as well as 
S = oo localized spins in infinite dimensions with use of 
the DMFT. It is found that quantum fluctuations of lo- 
calized spins partly destabilize the ferromagnetic ground 
state and the paramagnetic phase expands in the case of 
S = 1/2. The F state for 5 = 1/2 has rather small po- 
larization even for a strong coupling, such as J^S ^ 5.0. 
This is in contrast to the results for S = oo, where the 
ground states were almost fully polarized.^**' On the 
other hand the AF ground state is stabilized by quantum 
effects for weak coupling. A similar resiilt was recently 
reported for a model with degenerate orbitals.^^^ It is in- 
teresting to examine whether these results still hold for 
different choice of DOS (for non-interacting electrons) 
from the Gaussian one. 

More accurate calculations to determine the precise 
ground state phase diagram is left for a future study. 
The obtained phase diagrams for the 5=1/2 local- 
ized spin systems do not show good convergence at the 
boundary of the ferromagnetic phase. Furthermore the 
ground state was obtained under the assumption of the 
two-sublattice structure. We found that the phase tran- 
sition to an IC state occurs in the case of S = oo and we 
expect it to occur for 5 = 1/2 as well. 

It is interesting to study the case of 5 = 3/2 and com- 
pare the result with those for S = 1/2 and oo, since the 
magnitude of the localized spins are 3/2 in real mangan- 
ites. This is now under progress and will be reported 
elsewhere. 
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Fig. 1. The magnetic susceptibility x for Jy{S = 5.0 and n = 1.0. 
(a) Wave-vector dependence of x- (b) Temperature dependence 
of the antiferromagnetic susceptibility {X(q) = —1.0). ^(q) is 
defined by (3.1) in the text. 




Fig. 2. The susceptibility for JjjS' = 5.0 and n = 0.8. (a) Wave- 
vector dependence of the susceptibility, (b) Temperature depen- 
dence of the ferromagnetic susceptibility (X(q) = 1.0). 
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Fig. 3. Phase diagram for 5 = oo at /3 = 400. Antiferromagnetic, 
ferromagnetic and paramagnetic phases are denoted by AF, F 
and P, respectively. PS means the region where phase separation 
occurs. 
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Fig. 4. The magnetization for S = oo, f3 = 400 and n = 0.4 as a 
function of JhS. 
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Fig. 5. Wave-vector dependence of the spin susceptibility for 
JjjS = 1.0, n = 0.66 and (3 = 110. 




Fig. 6. Phase diagrams of the ground state derived by the ex- 
act diagonalization method with Afg = 6 (a) and TVg = 8 (b). 
AF, F and P denote antiferromagnetic, ferromagnetic and para- 
magnetic phases, respectively. PS means the region where phase 
separation occurs. 
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Fig. 7. The electron density n as a function of fj, in the ground 
state for JjjS = 5.0 and Afg = 8. Circles, squares and triangles 
indicate the antiferromagnetic, paramagnetic and ferromagnetic 
states, respectively. 




Fig. 8. The kinetic energy (a) and the magnetization (b) of the 
ground state as functions of Jfj5 for n = 0.4 and TVg = 8. Rhom- 
buses and squares indicate the paramagnetic and ferromagnetic 
states, respectively. 
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Fig. 9. Density of states for JjjS = 5.0 at n = 1.0 (a) and n 
0.78 (b), and for JjjS = 0.5 at n = 1.0 (c) and n = 0.73 (d). 
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